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We consider fermions on a 2D lattice interacting repulsively on the same site and attractively on 
the nearest neighbor sites. The model is relevant, for instance, to study the competition between 
antiferromagnetism and superconductivity in a Kondo lattice. We first solve the two-body problem 
to show that in the dilute and strong coupling limit the s-wave Bose condensed state is always the 
ground state. We then consider the many-body problem and treat it at mean-field level by solving 
exactly the usual gap equation. This guarantees that the superconducting wave-function correctly 
vanishes when the two fermions (with antiparallel spin) sit on the same site. This fact has important 
consequences on the superconducting state that are somewhat unusual. In particular this implies a 
radial node-line for the gap function. When a next neighbor hopping t' is present we find that the 
s-wave state may develop nodes on the Fermi surface. 
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I. INTRODUCTION 

Superconductivity is possible when carriers attract. 
The attraction in metals is normally due to the phonon 
coupling. The strong coulomb repulsion is normally well 
screened in metals and under some conditions it can be 
overcome by the attraction leading to superconductivity. 
A locally attractive (non-retarded) model, such as the 
BCS one, can thus give a good qualitative description of 
the superconducting state since the correlation length £ 
is much larger than both the screening length and the 
interparticle density. 

This is not the case when the carriers are non con- 
ventional. A relevant example is the single-band Kondo 
lattice with a RKKY interaction. For large Kondo cou- 
pling the ground state of this system is formed by local 
singlets: it is thus insulating at half-filling. At small 
doping the carriers are bachelor spins: they feel a hard 
core repulsion, since they cannot sit on the same site at 
the same time. The RKKY antiferromagnetic exchange 
interaction gives an attraction among holes in the sin- 
glet channel, and thus it could lead to a superconducting 
stateEm. But in this case the hard-core repulsion cannot 
be neglected in the description of the s-wave supercon- 
ductivity. More details on the Kondo lattice model are 
discussed in Reference ||. 

In this paper we consider a simple model where the role 
of the hard-core repulsion can be analyzed in a transpar- 
ent way. For simplicity and relevance to physical sys- 
tems we consider a two dimensional square lattice with 
hopping and attractive interaction among nearest neigh- 
bor sites. Different range for the interaction and the 
hopping produces interesting results, that will be ana- 



lyzed. Similar models have been studied by different 
authors. Ohkawa and Fukuyama proposed an ansatz 
to solve the mean-field equations for superconductivity 
with local repulsion in Kondo aystemsd. Later Micnas, 
Ranninger, and RobaszkiewiczQ and Bastide, Lacroix, 
and Rosa Simoestl, applied that technique to calculate 
the superconducting critical temperature. Similar equa- 
tions in a different context were considered by Aligia and 
coworkers B 



We will reconsider the problem at zero temper- 
ature and investigate the competition among differ- 
ent symmetries of the order parameter as one spans 
the crossovegpfrpajt-, BCS,— suoerconductivity to Bose 
condensationD&yyy'BBBtjy. In particular we will 
clarify an-apparent disagreement with the prediction of 
RanderiaLZI on the presence of bound states and super- 
conducting instabilities in 2D systems. We will also dis- 
cuss the unusual nature of the superconducting state: the 
main peculiarity is the presence of a radial line of nodes 
in k space for the gap function A(k). 

The plan of the paper is the following. In Section || we 
define in detail the model we are considering. In Section 
III we solve the two body problem and we discuss the 
relevance to the Bose limit. In Section |y| we set up the 
mean field equations and find the resulting phase diagram 
for the ground state. The competition between s- and 



d-wave superconductivity is analyzed in Section [y| In 
Section |v| the effect of a mismatch between the hopping 
and the interaction range is considered by introducing a 
hopping t' among next-neighbor sites. 



2 



II. THE MODEL 

We consider fermions hopping on a 2D square lattice 
of spacing a with nearest-neighbor sites hopping matrix 
elements t. (In some cases we will consider also next- 
neighbor hopping t' .) The fermions interact with a local 
repulsive interaction parameterized by U, that we will 
set to infinity. They furthermore feel an attraction V, 
which we assume to act only between nearest neighbors 
(in a real RKKY example that would not be the case) 
with opposite spins. The relevant parameters at zero 
temperature are only V/t and the filling of the band < 
n < 2. The spin does not play an essential role. Since 
we are interested in the singlet state, we can disregard 
the parallel-spin part of the interaction, the mean field 
equations would not be modified by that. The model 
Hamiltonian is thus: 

H = 5> k - M)4 ff c kCT + U ^(n lT - -){mi - \) 



-V^2(nn - h(n l+si - -) 



(1) 



where the vector S spans the four nearest neighbor sites 
of the square lattice, a is the spin projection in a chosen 
direction and t k = —2t[cos(k x a) + cos(fcj,a)] is the single- 
particle dispersion relation. 

The Hamiltonian in the form of Eq. (|l|) is particle-hole 
symmetric: for half- filling fi = 0. For convenience it can 
also be written in a non p-h symmetric, but shorter form: 



kcr i i,S 

(2) 

where \i Q = [i+(U— 4V)/2. The two forms are equivalent 
but from a constant. The interaction term can be written 
in Fourier space as usual: 



53 ^( k 2- k 3)c klT c L c : 

ki+k2=k 3 +k 4 



k 3 | C k 4 T 



(3) 



with V(k) = U - Vw(k) and w(k) = J^s e' H = 
2[cos(k x a) + cos(kya)]. 

In the following we investigate the superconducting or- 
der in the ground state of the Hamiltonian ([!]). In order 
to treat the strong repulsion correctly we consider first 
the two body problem, that is relevant in the dilute limit. 



choice between s- and d-wave is determined simply by 
the lowest 2-body bound state, if present. 

We thus begin our investigation by solving the two 
body problem: we find the threshold for the appearance 
of bound states. The two-body problem is exactly solv- 
able by diagonalization of a 5 by 5 matrix, corresponding 
to the central and 4 neighbor sites involved in the equa- 
tions of motions. This can be done directly in real space, 
but we will proceed in momentum space in order to use 
the same approach for the mean field equations in the 
next Section. In /c-space the simplification comes from 
the separability of the interaction potential. As a matter 
of fact, Vkk' = V(k — k'), defined above can be written 
in a separable form as follows: 



Vkk' = Y w a( k ) V "/3 w /3(k') 



(4) 



a/3 



where 



w (k) 
wi(k) 
u> 2 (k) 
w 3 (k) 
u>4 (k) 



= f 



[cos^a) + cos(fc a a)]/V2 
[cos^a) — cos^a)] /v2 
sin(fc 2 ,a) 
sm(k y a) 



(5) 



Vq/3 = v a 5 a p, with vq — U and v a = —2V for a ^ 
0. We have chosen a symmetrized and anti-symmetrized 
combination of cosines in order to exploit the symmetry 
of the lattice. 

The Green's function for the two-body problem on the 
basis |qQ) = c^ +Q/21 .cl q+ Q^JO) is diagonal in Q. We 

consider only the case Q = The corresponding equa- 
tion of motion is thus that of a particle in an external 
potential: 

(UJ - 2i k )£ k p(w) - 2J Vk; :q <? q pM = <5kp . (6) 



Since V is separable the above integral equation re- 
duces to an algebraic one for 5 parameters. The substi- 
tution of the separable form of V gives: 

(cj - 2i k )C/ kp (a;) - Y w a (k)A ap (uj) = 5 kp , (7) 



where 



A ap {u>) = v a Y w a( c L)Oc i p{oj) ■ 



(8) 



The Green's function is thus: 



III. TWO-BODY PROBLEM 

If the attractive interaction is strong enough to bind a 
pair we expect to reach the Bose condensed limit for a 
low density of carriers. It is thus crucial to know whether 
a bound state is present or not. In our case the hard-core 
repulsion will fight against an s-wave bound state, while 
it will not affect the <i-wave state. In the dilute limit the 



<W + Y Wa ( k )^"p M 



(9) 



where G±{w) = (u- 2i k )~ 1 . Substituting @ into (§) we 
obtain the following closure equation for A ap : 

A ap (uj) = v a w a (p)g°(u>) +v a y]B a p(u)App(ij) , (10) 
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FIG. 1: Two-body s-wave and d-wave bound state energy as 
a function of V/t and measured in terms of the band width 
2D — St. In both cases we find a threshold for the appearance 
of the bound state. 



where 



S a/9 (w) = 5^«; a (q3«; / ,(q)g°(w). (11) 



The set of linear equations (|T^) for A can be solved by 
inversion of the 5x5 matrix, 1 — V B(uj) 

A ap (co) = Y,T<*p{u)w p {v)G° p (w) (12) 

P 

where we have defined T(oj) = [1 — VS((jj)] _1 V. The final 
expression for the Green's function is: 



<5k P + g°(uj)^2T a p(uj)w a {k)wp(p) 



a/3 



(13) 

Note that T is the usual T-matrix of scattering theory: 



(14) 



a/3 



The problem can now be solved for any value of U, but 
since we are interested in the U = oo limit it is convenient 
to write T as: 



r = {v- 1 -b)- 1 



(15) 



where (V~ 1 )oo = 1/U = 0. In this way we eliminate U 
from the outset. One can readily verify that for U = oo 
the resulting local Green function Gfei = 0) actually 
vanishes. As it should, the hard-core repulsion prevents 
two fermions with opposite spin projection to sit on the 
same site. 



Let us come back to the explicit evaluation of T in 
two dimensions. The B matrix is in block form, C/° is in 

-k x and k y — > — k y . W2 is 



fact even under both k x 



odd under k x <— ► ky, and w 3 and W4 are odd under k x ► 
— k x and k y — > ~k y respectively. These symmetries leave 
few non vanishing elements and split the matrix in block 
form, corresponding to the irreducible representations of 
the symmetry group of the lattice: 



B 



/B 00 An \ 

B 01 B n 1 

B 22 

B 33 

V B 33 J 



(16) 



The bound states are given by the solution of the equa- 
tion det(V _1 — B) = 0. We obtain thus the following 
equations: 



— 1)1 R 
~ 7^ 1 



1 

2V 
1 

2V 



B01 
Boo 



— B21 = 



E 



[003(^0) — cos(k y a)] 2 



u> — 2ik 



1 

2V 



— = B. 



33 



E 



sin 2 (fca;a) 



(17) 
(18) 

(19) 



Eqs. @, @, and @ refer to s-, d-, and p-wave bound 
states, respectively. An additional s-wave bound state 
has been eliminated by setting U — 00. Note that U 
does not affect the p- and d-solutions, since the wave- 
function vanishes at the origin. In contrast, the s-wave 
bound state is strongly affected by U . As a matter of fact, 
Bqo(uj), -Boi(w), and Bu(w) all diverge logarithmically 
for uj — * 2D = 8t, with 2D the single particle band width. 
Using Eqs. (^0|) and (^l|) below one can verify that these 
divergences disappears in Eq. ([l7]), due to mutual can- 
cellation: the s-wave bound state is present only for V 
larger than a threshold V t s , while for a purely attractive 
potential binding always occurs in two-dimensions. 

It is not difficult to prove the following exact relations 
among the different £Ts by summing and subtracting uj 
and cosine terms in the numerator: 

ujB Q0 {uj) + 4V2 tB Ql (uj) = 1 , (20) 
ojB 01 (lu) +iV2tB n (uj) = 0. (21) 



Substituting Eqs. (gOj) and (|2l|) into Eq. 
bound state equation becomes: 



1l 71) the s-wave 



2t 
V 



UJ 

8t 



1 



8f B 00 {ujb) 



-wave , 



(22) 



where explicitly Bqo(uj) = X)k( w — 2ik) _1 - At the thresh- 
old V t s the bound state energy lies at the band bottom 
uj = —8t where -Boo diverges: hence the threshold V t s 
for a s-wave bound state is simply V = 2t. The s- 
wave solution appears non analytically at the threshold: 



ujb{V) 



-2D - 4iexp{-7rU t 7(T/ - V t s )}. For values 
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of V > V t s Eq. (|22j) can be solved numerically. We plot 
the result in Figure |l| together with the d-wave bound 
state energy obtained by numerical solution of Eq. (18). 
As it can be seen from the picture the <i-wave threshold 
V c d « 7.35 1 is higher than the s-wave one, (despite the 
fact that the latter is affected by the hardcore). This 
reflects the cost in kinetic energy (curvature of the wave 
functions) of changing the sign as we go around x = 0. 

We thus conclude that for V > V t s and for low den- 
sity the ground state of the system is a Bose condensate 
of fermions pairs in a s-wave bound state. The Bose 
limit is safely achieved only if the bound state energy 
\ll>b + 2D |, is much larger than the interparticle distance 
energy scale, i.e. ep ~ t n. Since ll>b + 2D vanishes ex- 
ponentially near V c the Bose region in the n-V diagram 
rapidly shrinks as V — > V c (cf. Fig. || in the following). 

This picture for dilute pairs is quite clear and it will 
not be modified by the introduction of a small t' (see 
following discussion). The opposite BCS limit is on the 
other hand more subtle. 



IV. MEAN FIELD EQUATIONS FOR 
SUPERCONDUCTIVITY 



In order to study the onset of superconductivity for 
small coupling and large density, we consider the usual 
BCS mean field theory. 

As in the two-body problem, the integral equation can 
be reduced to a simple algebraic equation. The s-wave 
solution is strongly modified by the hard-core interaction, 
and the order parameter has a k dependence that comes 
from the mixing of the wq and w\ terms, as in the two 
body problem. This mixing is crucial in order to ensure 
the vanishing of the superconducting wavefunction at the 
origin. 

The mean-field theory equations describing the super- 
conducting state of the Hamiltonian (|l|) are obtained by 
the usual decoupling: 



Hmf = &cL c k CT + [ AkC k T c -ki 



(23) 



The minimization procedure gives the familiar equations 
for the ground state: 



E 



in 



A k = -]Tv(k-k'; 



1Ev 



(24) 
(25) 



We concentrate now on the gap equation (|25[). We 
exploit again the separability of the interaction to write: 



A k r„ir, t (k) ^ w a (k') ^ k X 



k' 



= 22 A a w a (k) , 



where the new coefficients A satisfy: 



A 



Y,W a p[A}Ap, 



(26) 
(27) 



and we have defined the matrix W, analogous of the ma- 
trix B for the two-body problem: 



W a p[A] 



k 



w a (k)wp (k) 
2£ k 



(28) 



The A dependence comes through A in E. In matrix 
form Eq. (|^) reads simply (1 + VW)A = 0. As before 
it is convenient to write the equation for A in terms of 
V" 1 , in order to perform easily the limit U — ► oo. Eq. 
( ^7| ) thus becomes: 



(V" 1 + W[A])A = 0. 



(29) 



In order to find A k we need to solve the nonlinear 
equation (|2^), analogous of the linear equation (|To| ) for 
the two-body problem. Let Ai and v l a be the eigenvalues 
and eigenvectors of V -1 + W. We project A a —J2i c i y a 



on the v l a , so that 



(30) 



Cj is non vanishing only if Ai = 0. The corresponding 
coefficient c = Cj is fixed by the equation Aj[c] = 0. The 
form of the solution for A is simply A a = cv l a , and for 

A k = cE a "> a (k)- 

Due to the symmetry of the crystal, the matrix W 
is already in the form ( |30"| ) for a = 2, 3, and 4. The 
conditions = in these cases give the two equations 



1 

2V 

1 



[cos(fca;a) — cos(fcj,a) 



2D k 



— = y 

2V ^ 



sin 2 (k x a) 
2D k 



(31) 
(32) 



that correspond to the d- and p-wave solutions. In the 
"s" subspace {wq,wi} we are instead left with a 2 by 2 
matrix: 



where £ k = i k - /i + V(0)n/2, Eu 



and 



for consistency we retained the Hartree term. Putting 
together the Hartree shift and the original shift of the 
chemical potential we have £ k = t k — with // = fi + 
Tj(U — AV)(n — 1). (There is no Fock exchange term 
because of our interaction between opposite spins only.) 



v- 



W\A] = I Wo ° Wo1 
11 [ Wqi -1/2V + W n 



(33) 



The condition of vanishing of the determinant gives: 
I W$ 1 



2V Woo 



Wu, 



(34) 
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When a solution is found the form of the gap function 
will be given by the corresponding eigenvector of (33): 
v = (W i,-Woo), thus 



1 



cos fc x a 



cos k y a 



V2 



where 



r = W 00 /W 10 . 



(35) 



(36) 



The mixing of the two u>'s is crucial in order to take 
into account correctly the hard-core repulsion. This was 
already clear from the two-body problem, where the ex- 
act solution gives a wave function that is a superposition 
of wq and w\. At the MF level we could, in principle, 
search for solutions with only one of these two compo- 
nents different from zero: such a procedure is wrong. The 
simplest way to see that is to evaluate the anomalous cor- 
relation function JF(x — y) = (cj(x)cj(y)). This function 
is the equivalent of the wavefunction for the two-body 
problem: it gives the probability of finding two fermions 
on the same site, thus for infinite U it is bound to be zero 
at the origin. We can verify that our solution correctly 
gives F(0) = 0: 



^(«) = E^ = E 



2£ k 



Hi 



k K 



0(1 



E 



wi(k) 



(37) 



It is clear that if, out of the two order parameters, only 
one is present, .F(O) would be always non-vanishing. For 
p- and ei-wave, the symmetry alone is sufficient to guar- 
antee the vanishing of the wavefunction. 

The condition (|37| ) has important consequences on the 
superconducting order. It is clear that in order to ful- 
fill ( |37j ) Ak must change sign in the Brillouin zone. In 
particular, since wi(k) — — tk/(2-\/2i) 5 Ak is constant on 
surfaces of fixed energy: 



A k = A(a) = A[l + r(£k + n')/{2y/2t)] 



(38) 



Thus the line of nodes corresponds to a given energy £k 
instead of given direction, as it would be in p- or <i-wave 
superconductivity. The energy for which A vanishes is 
6v = -2y/2t/r - //. Other consequences of (|37| ) will 
be discussed in more details in Section IV D, after the 
solution of the self-consistent equations. 

Before proceeding, we find a simpler expression for the 
relevant Ws. Adding and subtracting p! in the numera- 
tor of W\§ we have: 



Win = - 



2V2t W °° 2x2/ 



—y 



2E k 



(39) 



We use Eq. 
We obtain: 



([24]) fixing n to eliminate the second term. 



2p'W 00 + 4V2tW ] 



10 



(40) 



Superconducting 




FIG. 2: Boundary between the s-wave superconductor and a 
normal phase. The dashed line is the line where A/(// +D) = 
.1, thus it gives an indication of the BCS region. The dotted 
line is the line where the chemical potential it exactly at the 
band edge: // = —D, and can be considered the beginning of 
the Bose region. The dot-dashed line in the dashed region is 
the analytical expression of Eq. ( fiij ) valid for small n. 



that is a generalization of Eq. (g0|) for the two-body prob- 
lem. Similarly for W\\ we have: 



8t 2 



W 0I 



M'(l - n) 



(41) 



where 



16t 2 4-^ #k 



16t 2 



16t 2 ^ £ k 

k 



In conclusion we obtain the following set of equations: 

(n-1) 2 



- = 26(//,A,r) 



16f 2 IFoo(M',A,r) 
n — 1 



2V2t AV2tW 00 (p',A, 



(43) 
(44) 



This form is more convenient for the following discussion 
since b remains finite when A — * 0. 

Equations (|3|) and together with Eq. (p4j), form 

the complete set of equations for the three unknowns 
{p',A,r} (the gap is defined in Eq. (p5|)). This will be 
our starting point to discuss the physics of the s-wave 
superconducting state. 



A. Threshold line for s-wave superconductivity 

We have seen that in the two-body problem a pair is 
bound only if V exceeds a threshold potential V t s . For 
small densities it is clear that for V > V t s the system is 
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superconducting and in the Bose limit, i.e. with a coher- 
ence length much smaller that the average interpacticle 
distance. As shown rigorously by Randeria et aZ.llH, in 
the low density limit and in 2D, the presence of a bound 
state in the two-body problem is a necessary and suffi- 
cient condition for s-wave superconducting instability to 
take place. We thus expect that for V < V t s and n — > 
the system is not superconducting. A separation line, 
starting at V t s for n — > 0, as function of n, separates the 
s-wave superconducting region and the normal behavior. 
The equation for this line of second-order transition can 
be found by setting A = in Eq. (|||) n being fixed by 
the chemical potential fj,' (—At < fi' < At = D). We 
obtain: 



1 

V 



8t 2 



1 



(45) 
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When A = 0, f/ + D = e° F is the Fermi energy of an 
ideal gas measured from the band edge. For small do ping 
e° F = 2ittn, and since for small n the integral in ( |4q ) 
is quadratic in ep we can neglect it and find a simple 
analytic expression for the threshold curve: 



= 2 



0(n 2 ). 



(46) 



The complete curve is shown in Fig. ||. The critical line 
increases smoothly from V/t = 2 for n = to V/t « 4.93 
at half- filling (n = 1). This can be easily understood 
by the fact that the particles feel more and more the 
presence of the hard-core repulsion. 



B. Near the threshold: BCS limit 

Near the transition line the system behaves like a BCS 
superconductor for weak coupling. In particular, the su- 
perconducting gap is much smaller than e F , the Fermi 
energy of the free gas, and /i' remains —D + e F . In order 
to prove this fact we assume that these properties hold, 
and we evaluate Eq. (Q) . Woo has the usual In A diver- 
gent behavior and it can be evaluated in leading order 
by assuming a constant density of states and constant 
Ak near the Fermi surface (tk = //). It is convenient to 
write the equation in terms of V — V t s (n) , by subtracting 
the critical line equation. In this way we obtain: 



A(/0 « 2/D 5 



/j,' 2 exp 



(1 - n) 2 VV t (n) 
m 2 [V -V t {n)] P (pi) 



(47) 

where /?(//) is the density of states per spin evaluated 
at //. Near the threshold A/e F <C 1 and thus the su- 
perconductor is of BCS type. Such a regime persists 
as long as A(//)/(// + D) is small. The dashed line 
in Fig. ^ indicates the values of n and V for which 
A(//)/(/F + D) = 0.1. The "BCS region" is restricted 
between that line and the critical line V t s - a narrow 
region compared to the remainder of the phase diagram. 
The behavior of the chemical potential is shown in Fig. ^. 



FIG. 3: Gap function at the Fermi energy A(/j,')/t as a func- 
tion of the density n, and for different values of the coupling 
V. Starting with the upper curve V/t is: 6, 5, 4.5, 4, 3.5, 3, 
and 2.5. One can easily see the square-root behavior for small 
n, typical of a Bose behavior. For V < 4.93, near the critical 
point value of the density a BCS exponential decay is evident. 
For V > 4.93 and at half-filling, the gap function A(fi') van- 
ishes on the Fermi surface remaining finite elsewhere. Since 
at half-filling the model fails, we do not discuss further this 
behavior. 



C. Bose limit and intermediate regime 

As discussed above, for V > V t s and small density we 
inevitably reach the Bose limit. This is confirmed by the 
mean field equations. Indeed for n — > the quantity 
A/l/ti' + D\ vanishes and // < —D. One can thus verify 
that the Eq. (^) reduces to Eq. fl22]) for the bound state 
energy of the two-body problem with \i' = ojb/2. At the 
same time A is fixed by Eq. (^) that gives A ~ tyfn. 
This behavior of A can be seen in Fig. || for small n. 
For intermediate couplings and densities the numerical 
solution crosses over between the BCS and Bose limits 
described above. In Fig. |^ we plot the lines of constant 
A(fi')/(fi' + D). We have seen that this ratio character- 
izes both the BCS (small positive value) and the Bose 
limit (small negative value). The phase diagram for the 
crossover from BCS to Bose condensation in the purely 
attractive case was given in Ref. [ji] for both s-wave (local 
attraction) and d- wave (nearest neighbour attraction) . A 
comparison of the results shows that the effect of the re- 
pulsion is to eliminate a large portion of the s-wavc BCS 
region, which becomes normal. The <i-wave case will be 
discussed in Section |v|. 

All these curves are plotted in the range of densities 
(0,1): it should be stressed that the have no real meaning 
near half filling since the mean field approach misses the 
insulating state. We expect that the approximation re- 
mains valid as far as ZkV, the weight of the quasiparticle 
peak, remains of the order of 1. We also remember that 
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FIG. 4: The chemical potential fi' measured from the bot- 
tom of the band —D, minus the free gas Fermi energy e° F . 
The curve are at fixed V/t which takes the values: 6, 5, 4.5, 
4, 3.5, 3, and 2.5 from the top to the bottom curve. The com- 
parison is particularly interesting near the critical value of n. 
There the behavior is BCS and the system is only slightly 
different from the free gas, since ef ~ (°f an d A < tf. For 
V larger that 4.93 no critical density exists, and the system 
never become strictly BCS. 




n 



FIG. 5: Phase diagram of the model for n and V/t, U = +oo. 
In the shaded region s-wave superconducting is not stable. In 
the rest of the diagram we find that s-wave superconductivity 
smoothly crosses from the BCS regime, near the critical line 
(the border of the shaded region), to the Bose regime, near 
the n = line. In order to characterize better this crossover 
we report the lines of fixed ratio A(/i')/(/x' + D). From the 
upper line to the critical line the ratio is cxj, 10, 5, 1, and 
.1. Surprisingly the BCS region is narrow, and the genuine 
crossover region covers most of the diagram. 



in the corresponding Kondo problem the limit n 
corresponds to no particles at all. 



D. Nodes of A(£) 

We have seen above that A depends only on £ and that 
it vanishes for £ = £jv- As a consequence we can write 
A(£) as follows: 



MO 



Ar 
2>/2t 



N 



(48) 



where, substituting Eq. ( p4| ) in the expression for £jv we 
have: 



1-n 



N 



2WooO',A,r) 



(49) 



For small value of A, £at ~ — t/ ln(A/i), thus it vanishes 
logarithmically in A. This means that approaching the 
threshold line V* the lines of nodes approach the Fermi 
surface. It turns out that keeping it at a logarithmically 
small distance (in A) suffices to fulfill the hard-core con- 
dition ©. 

The presence of a line of nodes near the Fermi sur- 
face modifies the usual shape of the occupation number 
distribution rik- Specifically, it introduces a particle-hole 
asymmetry for |£| w For large values of the energy 
the n distribution is anomalous. Usually, far from the 
Fermi surface, rik is either or 1, depending on the sign 



of £. In contrast our A(£) is ~ £ when |£| is large with 
respect to both A and £at, leading to 



»(0 - ~ 



1 - sgn(£) 



1 



y/l + A 2 r 2 /8i 2 



(50) 



This behavior can be recognized in Figure [3] where the 
form of n(£) is shown for V/t = 4.5 and for two different 
values of n. Inspection of that figure shows that 
remains symmetrical around the chemical potential only 
for |£| <C A(//). The asymmetry for |£| > £jv is clearly 
apparent in the inset. 

The aforementioned asymmetry explains the anoma- 
lous behavior of the chemical potential in the supercon- 
ducting state seen in Fig. ^. The chemical potential is 
normally shifted down with respect to the free gas value! 
The BCS distribution function is particle-hole symmet- 
ric, and states are "missing" below the lower band edge: 
/i must go down in order to retain the same density. In 
contrast here the lower spectral density near £ = £/v leads 
to an increase of /x in a narrow range where A and £yv 
are comparable, as shown in Figures || and [?]. 

V. d-WAVE SUPERCONDUCTIVITY 

We consider now the possibility of d-wave supercon- 
ductivity. By symmetry, the hard-core part of the inter- 
action is not seen by the two-body wave function, thus 
the scattering states should always feel an attraction at 
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FIG. 6: Occupation number n k = n(£ k ) for V/t — 4.5 and 
two different values of n: n = 0.4 (left curve) and n = 0.642 
(right curve). In the inset the region where n k and Ak van- 
ishes is shown enlarged. The dashed line indicates the value 
of the chemical potential, and the dotted one £at. 
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FIG. 7: The three relevant energy scales of the problem: 
(dashed lines) \fj,\, (dotted lines) £jv, and (full line )A(/z) are 
plotted for the same values of V of Fig.[| As expected, £n 
vanishes linearly in n — n t , while A vanishes exponentially in 
l/(n - n t ). 



finite energy. This is enough to induce BCS supercon- 
ductivity at any finite density and for any value of the 
coupling constant V. In contrast the Bose limit can be 
reached only if a bound state is present in the two body 
problem, V > V t — 7.35 1. Remember that in this limit 
s-wave always wins, we thus restrict to V < V t d . 

In the opposite limit, for small values of the coupling 
V <C V t s , the situation is again clear: we expect a stan- 
dard BCS d-wave superconductors, since there is no s- 
wave solution there. (The presence of a small hopping at 
nearest neighbors is discussed in the next section and it 
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FIG. 8: Ground state energy gain \5E\ for s-(full lines) and 
d-wave (dashed lines) superconducting state as a function of 
the density for different values of V (indicated in the figure) . 
The realized state is emphasized with a thicker line. 

will not change the final outcome.) 

The situation is much more interesting at intermediate 
couplings, V t s < V < V f d . In this region there is compe- 
tition between the two possible superconducting orders: 
the lowest in energy will be the actual ground state. We 
analyze qualitatively this competition as a function of the 
density n for a fixed value of the coupling V . For small 
density (n — > 0) the d-wave order parameter has a BCS- 
like form with an effective interaction V e ff = V (w^fs 
where (u>|) fs is the average over the Fermi surface of the 
angular factor. Since u>2(k) ~ (k* — ky), V e ff, which is 
of order w\, is reduced by a factor n 2 . The gap Ak is ex- 
ponentially small and the resulting energy gain vanishes 
extremely rapidly with n since SE « — (A^)irg/2. 

The s-wave solution, instead, is in the Bose limit, thus 
its energy gain is simply the number of bosons multiplied 
by the bound state energy: SE = naJs/2. For n — > the 
s-wave is always preferred to the d-wave. 

The opposite holds near the threshold for s- 
superconductivity: there d-wave still lowers energy, and 
is thus preferred. We expect a first order transition in be- 
tween. The question is "where"? Will the BCS s-region 
be completely sweeped out by the d-wave BCS state? 

In order to answer this question we need an explicit 
calculation of the ground state energy gain for the super- 
conducting state: 

5£ = 2^t k n k -2 *k-^A k JT k . (51) 

k k<k F k 

The last term in ( pl| ) is — A^/(2V^) for the d-wave sym- 
metry, and A 2 r 2 /(2V) for s-wave. Ad is obtained by 
solving Eq. Q and Eq. (||) with A k = A d w 2 (k). A 
and r are the quantities introduced in J35] ) and discussed 
in the previous Section. Numerical results are shown in 
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Fig. g. One can clearly see the crossing of the two en- 
ergy values indicating the first order transition. For small 
density, as expected, the s-wave state is always lower in 
energy (left side of the diagram) . Note however that the 
value of the energy for which the transition takes place 
depends strongly on the value of V considered. We can 
distinguish two cases: 

(i) V slightly above V t s (n = 0) the s-wave supercon- 
ductivity spans the whole crossover from the Bose limit 
(for n — ► 0) to the BCS one. The first order transition 
to the d-phase occurs deep in the s-BCS region. Note 
that this happens at tiny values of SE (for instance for 
V = 2.5 1, SE c /t < KT 10 ). Between the s- and the d- 
wave superconducting phases there exists a region that 
for all practical purposes is normal. Increasing the den- 
sity at fixed V one thus find that superconductivity is 
reentrant. 

(ii) At larger values of V, the s- and d-state energies 
cross for a value of the energy gain reasonably large (of 
the same order of the maximum value achieved along the 
s-wave evolution). Thus a nearly monotonic behavior of 
SE results (cf. the line V = 6i ). 

In conclusion we found that the d-wave competes with 
the s-wave only for large values of V. When V is near 
the threshold a quasi-normal region between the s-wave 
and the d-wave can survive, and the BCS-Bose crossover 
in s-wave can take place fully. 



VI. NEXT NEIGHBOR HOPPING AND d-WAVE 
s-WAVE COMPETITION 

In this Section we consider the effect of the inclusion 
of a nearest-neighbors hopping integral t' . The effect 
of this term on the superconducting critical tempera- 
ture has been investigated in Reffi and in more details 
in Ref.B. Here we consider ground state properties and 
in agreement with Refs.ou we find that the introduction 
of t' strongly modifies the s-wave phase diagram. As a 
matter of fact, the threshold line V t s is washed out by an 
arbitrary small t' . A transparent proof of that has been 
given in Ref.B. We reproduce here their arguments for 
zero temperature. Equation (|34| ) for A can be rewritten 
as follows: 

— \— = (w\) - ( Wl ) 2 = ((ioi - { Wl ) f) , (52) 

IV Woo 

where we have defined: 

k k 

When A is small, the weighted average over the Bril- 
louin zone is peaked on the Fermi surface. If w\ (k) has 
the same k dependence of the kinetic energy, the right- 
hand side of ( j52] ) vanishes at leading order in A, and the 
logarithmic divergence in Woo wm be ineffective: super- 
conductivity exists for V larger than a threshold Vt- In 



contrast, when the spectrum t^ is such that wiCk) varies 
on the Fermi surface, the right-hand side of (p2| ) is finite 
for A — > and superconductivity is possible however 
small V. n 

This is discussed extensively intl Here we want to 
point out two interrelated facts: (i) the relation of su- 
perconductivity with the presence of a bound state in 
the two-body problem, (ii) The unusual nature of the 
superconducting state created by a small t' in the region 
of the phase diagram where the system would be normal 
but for t' . |— - 

Randeria et alx-H proved, within a two-dimensional 
continuum model, that the existence of an s-wave bound 
state for the two-body problem is a necessary and suffi- 
cient condition for the existence of s-wave superconduc- 
tivity at low density. In the model at hand the presence 
of a small t' ^ makes superconductivity always possi- 
ble while it maintains a threshold for the two-body bound 
state. That threshold stems from the fact that w\ has no 
fluctuations near the band edge - hence no logws term 
in the pair propagatoc. This result seems to contradict 
the prediction of RefO. 

The solution of this inconsistency comes from the k 
dependence of A at the Fermi surface. In fact, even if 
the order parameter is of s-wave symmetry, it actually 
changes sign 8 times on the Fermi surface and it has a 
zero average on it. This can be verified by solving the set 
of Eqs. (|3l] ) and (|24) in power series of t' . One finds the 
following expression for A^ 

A k = A - t [„(k) - (v) FS ] + 0(t' 2 /^') (54) 
M 

where t^ = — 2t^/2wi(k) + t'v(k). More generally, it is 
clear that for any perturbation w(k) to the nearest neigh- 
bor hopping, Ak will change sign on the Fermi surface, 
and since it has s-wave symmetry this must happen at 
least 8 times. 

Due to the unusual behavior of Ak the results of Ran- 
deria et al. do not apply to this s-wave state which is, 
in practice, close to an angular momentum I = 4 state 
(which is isotropic as regards cubic symmetry.) 

Hence s-wave superconductivity survives the hard core 
for arbitrary small values of the attraction. But since 
the superconducting state is close to I = 4, the resulting 
value of A (or T c ) is extremely small, even compared 
to the e?-wave solution. We checked explicitly this fact 
for our model. We considered values of t' < t/2, since 
at this value the Fermi surface becomes double sheeted. 
We calculated the effective s-wave coupling 

V eff/ V = (»?) " M 2 - [(v 2 ) ~ (v) 2 ] t' 2 /(8t*) (55) 

as a function of the Fermi energy and of the coupling t' . 
We compared Vj!ff with the effective ci-wave coupling: 
V*ff/V = (w%). In all cases we found that V&f > Veff- 
We thus conclude that the introduction of a t' < t/2, 
or of other similar perturbations, does not change quali- 
tatively the above discussion, since d-wave superconduc- 
tivity always hides this anomalous s-wave state. 
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FIG. 9: Phase diagram considering the competition between 
s- and d-wave. The thick line corresponds to the first-order 
phase transition between the s- and d-wave superconducting 
order. On the s-side we indicate the BCS and Bose conden- 
sation regions. On the d-side the region where SE/t < 10~ 6 
is labeled as "quasi-normal". 



VII. CONCLUSIONS 

We have studied the appearance of superconductivity 
on a 2D lattice in presence of a hard-core repulsion and of 
a nearest neighbor attraction. We constructed the mean 
field solution and compared it to the two-body problem in 
vacuum. Even if the procedure is restricted to mean field, 
it correctly forbids the double occupancy of the same site 
of the superconducting wavefunction. The main results 
are the following: (i) the s-wave solution is suppressed 



for small values of the coupling at any density. The in- 
troduction of additional small hopping integrals cannot 
change this, since the <i-wave solution is always preferred, 
(ii) the d-wave solution is always possible, but the actual 
value of the energy gain becomes extremely small at low 
density. Thus for small V and n the critical tempera- 
tures are tiny: in practice at any temperature the system 
is normal, (iii) For V larger than the threshold for the 
s two body-bound state the system exhibits a crossover 
from the Bose-Einstein condensation of fermions pairs to 
a BCS behavior as the density is increased. A first order 
transition to the <i-wave state occurs. The situation is 
summarized in the phase diagram of Fig. |^. 

When only nearest neighbor attraction is present the 
competition between e?-wave and s-wave at T c Jcms been 
considered recently by Wallington and AnnettE3. They 
found that for any coupling the s-wave has the lowest 
critical temperature, but for a small region near half- 
filling, where the van Hove singularity of density of states 
weighted with the angular factors stabilizes the d-wave 
phase. The presence of the hard-core repulsion allows d- 
wave superconductivity to appear on a much larger por- 
tion of the phase diagram. The correct treatment of the 
repulsion is thus crucial to study the competition between 
s- and d-wave. 

The technical procedure developed can be applied to 
the Kondo-lattice problem sketched in the introduction. 
However, one has to keep in mind that in that prob- 
lem superconductivity and Kondo hybridization should 
be considered on the same footing, since each of the 
two order parameter reacts on the other. We will 
not try thus to apply the above result directly to the 
Kondo lattice model: we leave the detailed solution of 
the superconductivity-Kondo hybridizatipa in presence 
of the hard-core repulsion for future workEl 
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